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Janus particles have attracted significant interest as building blocks for complex materials in recent 
years. Furthermore, capillary interactions have been identified as a promising tool for directed 
self-assembly of particles at fluid-fluid interfaces. In this paper, we develop theoretical models 
describing the behaviour of magnetic Janus particles adsorbed at fluid-fluid interfaces interacting 
with an external magnetic field. Using numerical simulations, we test the models predictions and 
show that the magnetic Janus particles deform the interface in a dipolar manner. We suggest how to 
utilise the resulting dipolar capillary interactions to assemble particles at a fluid-fluid interface, and 
further demonstrate that the strength of these interactions can be tuned by altering the external 
field strength, opening up the possibility to create novel, reconfigurable materials. 

PACS numbers: 47.ll.-j, 47.55.Kf, 77.84.Nh. 


I. INTRODUCTION 

Golloidal Janus particles have drawn special attention 
during the past two decades for their potential in ma¬ 
terials science [T]. Janus particles are characterized by 
anisotropic surface chemical (e.g. wetting or catalytic) 
or physical (e.g. optical, electric, or magnetic) proper¬ 
ties at well-defined areas on the particle. This combi¬ 
nation of chemical anisotropy and response to external 
fields makes Janus particles promising building blocks of 
reconfigurable and programmable self-assembled struc¬ 
tures [2HS]- 

Janus particles strongly adsorb at fluid-fluid inter¬ 
faces [7], making the formation of 2-D structures ac¬ 
cessible. For symmetric Janus particles composed of 
hydrophobic and hydrophilic hemispheres, the equilib¬ 
rium contact angle is 90° since each hemisphere im¬ 
merses in its favourable fluid, and the interface re¬ 
mains flat [5]. However, due to surface roughness [S], 
anisotropic shape nni im, or the influence of external 
forces, Janus particles can tilt with respect to the in¬ 
terface. In a tilted orientation, the fluid-fluid interface 
around the Janus particle deforms in a dipolar fash¬ 
ion in order to fulfil boundary conditions stipulated by 
Young’s equation [12]. Assuming small interface de¬ 
formations, the particle-induced interface deformations 
obey = 0, where h is the interface height, which 
can be solved using a multipolar analysis, analagous to 
2D electrostatics [T3|. These particle induced deforma- 


* q.xiel@tue.nl 
1 g.davies.ll@ucl.ac.uk 
t f.s.guenther@tue.nl 
S j.harting@tue.nl 


tions, called capillary interactions, can cause particles to 
attract and repel in specific orientations, making them 
a useful tool for controlling the behaviour of particles at 
interfaces. 

Previous investigations into capillary interactions be¬ 
tween particles at fluid-fluid interfaces have focussed 
mainly around two themes: particle weight-induced de¬ 
formations, which lead to monopolar interactions be¬ 
tween particles and are responsible for e.g. the Cheerios 
effect |14j ; and surface roughness or shape anisotropy in¬ 
duced deformations, which lead to quadrupolar interac¬ 
tions between particles [13] and are responsible for e.g. 
the suppression of the coffee ring effect [TB]. With re¬ 
spect to Janus particles, Brugarolas et al. m showed 
that quadrupolar capillary interactions induced by sur¬ 
face roughness can be used to form fractal-like structures 
of Janus nanoparticle-shelled bubbles. However, a signif¬ 
icant limitation of the above mentioned capillary inter¬ 
actions is that they are not dynamically tunable because 
they depend on the particle properties alone. 

Davies et al. m recently found a way of creating dy¬ 
namically tunable dipolar capillary interactions between 
magnetic ellipsoidal particles adsorbed at an interface 
under the influence of an external magnetic field. The 
structures that form depend on the dipole-field coupling, 
which can be controlled dynamically m- However, it 
is also desirable to create tunable capillary interactions 
between spherical particles without relying on particle 
shape anisotropy. 

In this paper, we show how to create tunable dipolar 
capillary interactions using spherical Janus particles ad¬ 
sorbed at fluid-fluid interfaces. The Janus particles have 
a dipole moment orthogonal to their Janus boundary and 
are influenced by an external magnetic field directed par¬ 
allel to the interface. The held causes the particles to ex- 
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perience a magnetic torque, but surface tension opposes 
this torque, and the particles therefore tilt with respect 
to the interface. When tilted, the particles deform the 
interface in a dipolar fashion. 

We develop a free energy based model of the behaviour 
of a single particle at the interface, and a model which 
takes into account small interface deformations. We nu¬ 
merically investigate the predictions of the models us¬ 
ing lattice Boltzmann simulations and highlight that the 
dipolar interface deformations lead to novel and inter¬ 
esting particle behaviour at the interface. Finally, we 
explain how to use these interface deformations to manu¬ 
facture tunable capillary interactions between many par¬ 
ticles at a fluid-fluid interface, which could be utilised to 
assemble reconfigurable materials. 

This paper is organised as follows. We present our 
hybrid molecular dynamics-lattice Boltzmann simulation 
method in Section |ll] In Section [TiTl we develop two the¬ 
oretical models describing the behaviour of Janus parti¬ 
cles at fluid-fluid interfaces. Section Hvl contains our sim¬ 
ulation results, and we compare these results with our 
theoretical models from Section [ITT] Finally, Section |V] 
concludes the article. 


II. SIMULATION METHOD 
A. The multicomponent lattice Boltzmann method 

We use the lattice Boltzmann method (LBM) to simu¬ 
late the motion of each fluid. The LBM is a local meso¬ 
scopic algorithm, allowing for efficient parallel implemen¬ 
tations, and has demonstrated itself as a powerful tool 
for numerical simulations of fluid flows ED]. It has been 
extended to allow the simulation of, for example, multi¬ 
phase/multicomponent fluids EUEH and suspensions of 
particles of arbitrary shape and wettability EDHUj. 

We implement the pseudopotential multicomponent 
LBM method of Shan and Chen |2T| with a D3Q19 lat¬ 
tice ED] and review some relevant details in the follow¬ 
ing. For a detailed description of the method, we refer 
the reader to the relevant literature. EH lUl I27ti2^ Two 
fluid components are modelled by following the evolution 
of each distribution function discretized in space and time 
according to the lattice Boltzmann equation: 

/^(f -b c* At, t + At) = f^{x, t)-^ [fiix, t) 

( 1 ) 

where i = 1,..., 19, /f(x, t) are the single-particle dis¬ 
tribution functions for fluid component c = 1 or 2, 
Ci is the discrete velocity in tth direction, and t‘^ is 
the relaxation time for component c. The macro¬ 
scopic densities and velocities are defined as p‘^{x,t) = 
Po^ifi{x,t), where po is a reference density, and 
u^{x,t) = i), respectively. Here, 

f^'^{p‘^{x,t),u^{x,t)) is a third-order equilibrium distri¬ 
bution function. When sufficient lattice symmetry is 


guaranteed, the Navier-Stokes equations can be recovered 
from Eq. Q on appropriate length and time scales |20j . 
For convenience we choose the lattice constant Ax, the 
timestep At, the unit mass po and the relaxation time r'^ 
to be unity, which leads to a kinematic viscosity = | 
in lattice units. 

The Shan-Chen multicomponent model introduces a 
mean-field interaction force 

Fc(x,t) = -^'''(x,t)^5cc' (x ,t)(x-x)(2) 

c' x' 

between fluid components c and c'm, in which x' denote 
the nearest neighbours of lattice site x and pcd is a cou¬ 
pling constant determining the surface tension. 'L‘^(x, t) 
is an “effective mass”, chosen with the following func¬ 
tional form: 

«'"(x, t) = «'(p"(x, t)) = l- (3) 

This force is then applied to the component c by adding 
a shift Au‘^(x, t) = to the velocity u^{x, t) in the 

equilibrium distribution. The Shan-Chen LB method is 
a diffuse interface method, resulting in an interface width 
of « 5Ax El]- 

B. The colloidal particle 

The trajectory of the colloidal Janus particle is up¬ 
dated using a leap-frog integrator. The particle is dis¬ 
cretized on the fluid lattice and coupled to the fluid 
species by means of a modified bounce-back boundary 
condition as pioneered by Ladd and Aidun ESI ISO] • 

The outer shell of the particle is filled with a “virtual” 
fluid with the density 

= p^{x,t)+ \Ap\, (4) 

Pvirt(^>0 = f{x,t)-\Ap\, (5) 

where p^{x,t) and p'^{x,t) are the average of the den¬ 
sity of neighbouring fluid nodes for component 1 and 2, 
respectively. The parameter Ap is called the “particle 
colour” and dictates the contact angle of the particle. A 
particle colour Ap — 0 corresponds to a contact angle of 
dp = 90°, i.e. a neutrally wetting particle. In order to 
simulate a Janus particle, we set different particle colours 
in well defined surface areas corresponding to the differ¬ 
ent hemispheres of the particle. 

We choose a system size S = 128 x 64 x 128 to eliminate 
finite size effects. We fill one half of the system with 
fluid 1 and the other half with fluid 2 of equal density 
(pi = P 2 = 0.7) such that a fluid-fluid interface forms 
at 2 / = 32. The interaction strength in Eq. is chosen 
to be gi 2 = 0.1 and the particle with radius i? = 10 is 
placed at the interface. We impose walls with mid-grid 
bounce back boundary conditions at the top and bottom 
of the system parallel to the interface, while all other 
boundaries are periodic. 
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(a) Equilibrium orientation (b) Tilted orientation 

FIG. 1: A single Janus particle adsorbed at a fluid-fluid interface in its equilibrium orientation (a) and in a tilted 
orientation (&). The Janus particle consists of an apolar and a polar hemisphere. The particle’s magnetic dipole 
moment m Is orthogonal to the Janus boundary, and the external magnetic field, H, is directed parallel to the 
Interface. The tilt angle is defined as the angle between the particle’s dipole moment and the undeformed 
Interface normal. Aa 2 and Api are the surface areas of the apolar hemisphere immersed in fluid 2 and the polar 
hemisphere immersed in fluid 1, respectively. The bold red line represents the deformed Interface and Ci is the 

maximal Interface height at the contact line. 


III. THEORETICAL RESULTS 

We consider a spherical Janus particle composed of ap¬ 
olar and polar hemispheres adsorbed at a fluid-fluid inter¬ 
face, as illustrated in Fig. The two hemispheres have 
opposite wettability, represented by the three-phase con¬ 
tact angles 9a = 90° -I- /? and Op = 90° — /3, respectively, 
where /3 represents the amphlphillcity of the particle. A 
larger /? value corresponds to a greater degree of par¬ 
ticle amphiphlllclty. In Its equilibrium state, the Janus 
particle takes an upright orientation with respect to the 
interface with its two hemispheres totally Immersed In 
their favourable phases, as shown in Fig. [la| 

The free energy of the particle In Its equilibrium con¬ 
figuration is 

= 7i 2A“ + 7aiA“‘ + 7 p2A“‘, (6) 

where "fij are the interface tensions between phases i and 
j and Aij are the contact surface areas between phases 
i and j, where i,j = {1: fluid, 2: fluid, a: apolar, p: 
polar}. For a symmetric amphiphilic spherical particle, 
the apolar and polar surface areas are equal Aai = Ap 2 = 

After switching on the horizontal magnetic field, H, 
the particle experiences a torque r = m x H that at¬ 
tempts to align the particle dipole axis with the field. 
However, surface tension resists the rotation causing 
the particle to tilt with respect to the Interface for a 
given dipole-field strength B = |m||H|, as Illustrated In 


Fig. Ib The tilt angle p Is defined as the angle be¬ 
tween particle dipole-moment and the undeformed In¬ 
terface normal (l.e. the y-axls). The Interface deforms 
around the particle so that each fluid contacts a larger 
area of its favourable particle surface (Fig. |Ib[ ). This In¬ 
terface deformation increases the fluid-fluid interface area 
and decreases the surface of each hemisphere contacting 
its unfavourable fluid. The free energy of the system 
Is reduced In total due to the dominant contribution of 
particle-fluid interface energies [T^] . The free energy of a 
Janus particle in a tilted orientation can be written as 


= 712A«“ + 7aiA«}‘ + 7p2A}i^‘ 

+ 7 a 2 Aa 2 + IpiApi J- B sln p. 


( 7 ) 


The free energy difference between the tilted orientation 
state and the Initial state is 


AA — Atilt ~ Aint 

= 712 - A“) + 7ai (A“* - A“‘) 

+7p2 (^p2* ~ ^p2 ) + 702^02 -b 7pl^pl 
+Bsm(p. ( 8 ) 

Since -b and A“* = A*!}* -b A^i, we obtain 

AA = 712 AA 12 -b (7a2 — 7al)Aa2 

+ (7pi ~ 7 ^ 2 ) Api -b Bsinp, (9) 

where AA 12 = A}'}* — A'^ is the increased fluid-fluid in¬ 
terface area. The particle obeys Young’s boundary con- 
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(a) (b) 


FIG. 2: Snapshots of a Janus particle at a fluid-fluid interface at (a) initial equilibrium state {ip = 0°) and (&) tilted 
orientation state {ip = 90°) as obtained from our simulations. The three-phase contact line undulates around the 

tilted particle so that the interface is deformed. 


ditions [8] 

COsOa = -, COS Op = — -(10) 

7l2 7l2 

For two hemispheres with opposite wettabilities, we ob¬ 
tain cos9A = —cosOp = — sin/3. Eq. (|^ is further sim¬ 
plified to 

AE = 7 i 2 AAi 2 -l- 7 i 2 (^a 2 + ^pi) sin/3-|-i?siny).(ll) 

Under the assumption of a flat interface [TU], ^^12 = 0, 
2^02 = ^pi = ^47ri?^ = 2ipR^. Therefore, Eq. ^ finally 
reduces to 


AE = 4(pi?^7i2 sin j3 B sin ip. 


( 12 ) 


There is no exact analytical expression for the free en¬ 
ergy of a tilted Janus particle at an interface that includes 
interface deformations, due to the difficulty in modelling 
the shape of the interface and position of the contact 
line. However, in the limit of small interface deforma¬ 
tions |13) , we derive such an analytical expression for the 
free energy. 

We consider micron-sized particles with a radius much 
smaller than the capillary length such that we can ne¬ 
glect the effect of gravity. We assume the pressure drop 
across the interface to be zero, leading to vanishing mean 
curvature according to the Young-Laplace equation. The 
mean curvature can be approximately written in cylindri¬ 
cal coordinates as m 


where h is the height of the interface. The radial dis¬ 
tance r and the polar angle O are defined with respect 
to a particle centred reference frame. Using a multipole 
analysis m, the solution of Eq. yields 


= Rm{r)^m{'&) (14) 

m>0 

oo 

^ \m 

Cm cos {m{'0 - -dm. o)) , (15) 

m>0 


where is the maximal height of the contact line, i?rn,o 
is the phase angle and is the radius where the particle 
and fluid interface intersect. Cc is approximately the par¬ 
ticle radius in the limit of small interface deformations. 
The monopolar term m = 0 is ommitted because we fo¬ 
cus on micron-sized particles where gravitational effects 
can be neglected. The dipolar term m = 1, results from 
an external torque on the particle, causing symmetric in¬ 
terface rise and depression around it. m = 2 denotes the 
quadrupole term, which dominates in the absence of any 
external forces or torques on the particle. We consider 
only the leading order m = 1 dipole term 

/i(r,i?) = Cl cos (i9 - di,o)—, (16) 

r 

where Ci is the maximal height of the contact line. We 
calculate the increased fluid-fluid interface area AH 12 by 
considering an infinitesimal element dA* = dx* x dy* 
of the deformed interface. In a local coordinate system 
rotated such that the slope is maximized along the y 
coordinate, we have [13] 


dx* = dx, 


dy* = y'dy'^ + dh'^ dy (^1 + ^(V/i)^^ , 


resulting in 


d{AAi 2 ) = dx* X dy* — dx x dy 

1 9 

= dxdy- (Vft.) 


or 


-I pOO p2T^ 

AH12 = - / / {'^hfrdMr. 

2 J r—Tc J '&—0 

In the dipole approximation, one has 

(Vk)-(V/.) = 

= drlr-^ 


(17) 

(18) 


(19) 


( 20 ) 


( 21 ) 







5 


and we obtain 

( 22 ) 

The areas Aa 2 and Api can then be written as 

pTT 

Aa 2 = Api = 2R^ip - h{r = reced'd 
do 

= 2i?V-2rcCi. (23) 

We assume r,. = R and the free energy taking into ac¬ 
count small interface deformations can then be written 
as 

A-E = ^ 7 i 2 Ci+ 4(i?^‘/5-i?Ci)7i2 sin/3-|-i?sin(/3.(24) 

We discuss the predictions of our models and compare 
them with our simulation results in the following section. 


IV. SIMULATION RESULTS AND 
COMPARISON TO THEORY 



FIG. 3: Reduced torque t jAp^x^ as a function of tilt 
angle for (3 = 21° (circles) , /3 = 30° (triangles) and 
P = 39° (diamonds), where Ap = ttR^. The symbols are 
simulation data and the solid lines represent hyperbolic 
tangent functions to fit the data. The fitted functions 
are integrated in order to obtain the free energy. For all 
amphiphilicities, the torque is linear for small rotations 
around the equilibrium position, before reaching a 
constant value at very large tilt angles ip —)■ 90°. 


In addition to the theoretical models we developed in 
Section III we numerically calculate the free energy of 


a single Janus particle adsorbed at a fluid-fluid inter¬ 
face using lattice Boltzmann simulations. Our lattice 
Boltzmann simulations are capable of capturing inter¬ 
face deformations fully without making any assumptions 
about the magnitude of the deformations or stipulating 
any particle-fluid boundary conditions. 

Fig. shows the initial equilibrium configuration of 
our Janus particle simulation, where each hemisphere to¬ 
tally immerses in its corresponding favourable liquid so 
that the interface remains flat. Fig. shows how the 
three-phase contact line and interface deform around the 
Janus particle as it tilts with respect to the interface. 

In order to obtain the free energy of the Janus particle 
as a function of tilt angle from our simulations, the total 
surface area of the deformed fluid-fluid interface and the 
corresponding particle surfaces have to be measured m 
1181131] . In this paper, we employ a simple method, which 
utilises the ability to easily measure the force applied to 
the particle by the fluids using the lattice Boltzmann 
method. 

We first determine the contribution to the free energy 
neglecting the dipole-field contributions by integrating 
the torque on the particle as the particle rotates quasi- 
statically, 


AE = / (25) 

^0 

where ptut is the tilt angle of interest. To do this, we 
rotate the particle on the interface until it reaches the 
desired tilt angle and then fix the position of the particle 
and let the system equilibrate. The remaining torque on 
the particle is the resistive torque applied to the particle 
from the fluid-fluid interface. 

Fig. shows the evolution of this torque versus 
the tilt angle for different amphiphilicities /3 = 21° (cir¬ 
cles), P = 30° (triangles), and P = 39° (diamonds), cor¬ 
responding to particle colours Ap = 0.10, 0.15, 0.20, 
respectively. For all amphiphilicities, the torque in¬ 
creases linearly as the particle rotates for small tilt an¬ 
gles, ip < 30°. As the tilt angle increases further p —> 90° 
the torque tends to a nearly constant value. We fit the 
torque with a hyperbolic tangent function, and inte¬ 
grate the fitted function numerically to obtain the free 
energy. 

In order to calculate the free energy of our small in¬ 
terface deformation based model (Eq. ([^), we measure 
the corresponding maximal height of the contact line 
as a function of tilt angle, as shown in Fig. Similarly 
to Fig. the height of the contact line increases lin¬ 
early for small tilt angles, and then reaches a plateau for 
large tilt angles. The height plateau demonstrates that 
the deformed interface area remains constant at large tilt 
angles. This indicates that at large tilt angles, only the 
particle-fluid surface energy, which increases linearly with 
increasing tilt angle, contributes to the change in free en¬ 
ergy, and therefore the torque r oc d{AE)/ dp is constant 
at large tilt angles, in agreement with Fig. 
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FIG. 4: Reduced maximal height of the deformed 
interface Ci/R as a function of tilt angle (p for 
/3 = 21° (circles) , /3 = 30° (triangles) and j5 = 39° 
(diamonds). The data (symbols) are fitted with a 
hyperbolic tangent function (solid lines). The maximal 
interface height increases linearly with tilt angle for 
small angles, and becomes constant for large angles. 


In Fig. we compare the free energy models that we 
developed in Section |III| assuming no interface deforma¬ 
tions (Eq. (12), dashed lines) and small interface defor¬ 
mations (Eq. (24), solid lines) with our lattice Boltzmann 
simulation data (symbols), which incorporates interface 
deformation fully, by measuring the free energies as a 
function of particle tilt angle as described above. 

Our undeformed interface model (Eq. (EH)) predicts 
that the energy varies linearly with the particle tilt angle, 
p, for all particle amphiphilicities, /3 = 21°, 30°, and 39°. 
The model also predicts that the free energy increases as 
the amphiphilicity increases from /3 = 21°to/3 = 39° for 
any given tilt angle. Our analytical model assuming small 
interface deformations (Eq. @) shows some interesting 
qualitative behaviour. Firstly, for small tilt angles, p < 
30°, the energy varies quadratically with the tilt angle 
for all amphiphilicities. This is because for small tilt 
angles, the maximal interface height varies linearly 
with the tilt angle (Fig. |^, ()i « Rp, and the free energy 
in Eq. (24) becomes approximately AE = 

Secondly, for larger tilt angles, p > 45°, the energy 
varies linearly with the tilt angle, due to the fact that for 
large angles, Ci is constant (Fig. |^. Therefore, Eq. (24) 
becomes AE = 4i?^7i2 sin Pp + C, where C is a constant, 
which explains the linear behaviour of the free energy for 
large tilt angles. 


FIG. 5: Free energy as a function of tilt angle p for 
P = 21° (circles), P = 30° (triangles) and P = 39° 
(diamonds), which are calculated using the analytical 
model that excludes interface deformations (Eq. (12), 
dashed lines), our analytical model which assumes small 
interface deformations (Eq. ( [24| , solid lines) and 
numerically integrated results (symbols) without 
considering a magnetic dipole-held contribution 
{B = 0). The undeformed interface model Eq. (12) 
shows non-negligible deviation from the simulation 
data. For small P and small p, the small deformation 
model Eq. (24) is in good agreement with the 
simulation results. 


The above analysis also explains why, for tilt angles 
p < 30°, the energy only weakly depends on the am¬ 
phiphilicity P of the particles: the amphiphilicity term 
sin P is only signihcant for large tilt angles. 

When comparing these two models with our simulation 
data (symbols), we hnd that the small deformation model 
captures the qualitative features of the data extremely 
well. In addition, the model quantitatively agrees with 
the numerical results for small tilt angles p < 30° for 
all amphiphilicities. As the tilt angle of the particles 
increases the quantitative deviation between the model 
and the data becomes more significant. However, for 
small amphiphilicities /3 = 21° the model and numerical 
data are in good agreement. 

In contrast, our theoretical model that takes into ac¬ 
count only the free energy differences between the par¬ 
ticle as a function of its orientation and therefore ne¬ 
glects interface deformations (Eq. El) performs much 
worse. The model captures the qualitative linear be¬ 
haviour of the numerical data only for large tilt angles 
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dipole-field strength B/ApY-|2 


FIG. 6 : Tilt angle tp as a function of dipole-field 
strength for different amphiphilicities /3 = 21°(circles), 
/3 = 30° (triangles) and j3 = 39° (diamonds). We 
compare the simulation data (symbols) with the tilted 
angles predicted using the analytical model which 
excludes interface deformations (Eq. (12), dashed lines), 
and our analytical model which assumes small interface 
deformations (Eq. (24), solid lines). Eor small 
dipole-field strengths, the tilt angle predicted by the 
undeformed interface model shows large deviations from 
the simulation data, whereas the tilt angle obtained 
from the small deformation model agrees well with the 
simulations. 


(/? > 50°, but with large quantitative differences that in¬ 
crease as the particle amphiphilicity increases. The re¬ 
sults in Eig. clearly show that interface deformations 
strongly affect the behaviour of a tilted Janus particle 
adsorbed at a fluid-fluid interface, and we note that our 
analytical model that includes small interface deforma¬ 
tions is clearly able to capture this qualitative behaviour. 

To verify the predictions of the analytical models in¬ 
cluding dipole-field contributions, we switch on the mag¬ 
netic field and numerically measure the time-averaged tilt 
angle of the particle after the system has equilibrated, as 
per Eq. (25). We obtain the tilt angle as a function of 
the dipole-field strength by minimizing the free energies 
in Eq. ([T^ and Eq. (241 with respect to the tilt angle. 


In Fig. we compare the predicted tilt angles from 
the free energy models assuming no interface deforma¬ 
tion (Eq. ( [T^ , dashed lines) and small interface deforma¬ 
tion (Eq. ( |24[ ), solid lines) with our numerical simulation 
data (symbols), which incorporates interface deformation 
fully, by measuring the tilt angle as a function of dipole- 
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FIG. 7: Plot of the relative height of the interface. 

Influenced by a magnetic torque, the Janus particle 
reorients, leading to dipolar interface deformations: the 
interface is depressed on one side of particle (black) and 
raised on the other side (yellow). The interface is flat in 
the orange/red regions. These dipolar interface 
deformations are dynamically tunable, providing a 
route to generate tunable capillary-driven assembly at 
fluid interfaces. 


field strength, for (3 = 21°, 30°, and f3 = 39°. 

For large dipole-field strengths B > 1 . 5 Ap 7 i 2 , both 
analytical models perform well by predicting the numer¬ 
ically measured tilt angles for all particle amphiphilic¬ 
ities p. This is due to the fact that large dipole-field 
strengths cause large tilt angles at which the deformed 
interface area stops increasing, as shown in Fig. In 
this regime, the equilibrium orientation only depends on 
the particle properties (radius R and amphiphilicity /3), 
interface tension 712 , and dipole-field strength B, which 
are incorporated into both analytical models. 

For small dipole-field strengths, the tilt angle predicted 
by the undeformed interface model Eq. (12) shows large 
deviations from the simulation data for all amphiphilic¬ 
ities. The deviations increase with increasing particle 
amphiphilicity. In addition, the undeformed interface 
model predicts a large critical dipole-field strength at 
which the particle begins to rotate, which we do not 
observe in the simulations. This critical dipole-field 
strength increases with increasing particle amphiphilic¬ 
ity. From Eq. (12), we can obtain the torque, r oc 
= 477^712 sin P + B cos ip where the first term is 


independent of the tilt angle p. Therefore, in the zero 
dipole-field state p = Q°, a, critical dipole-field strength 
Be = 477^712 sin P is needed to overcome this constant 
resistive component of the torque. 














Our small deformation model (Eq. ( [2^ ) shows signif¬ 
icant improvement at predicting the tilt angle compared 
with the undeformed interface model. The predicted tilt 
angle from this model shows good agreement with simu¬ 
lation data for weak dipole-field strengths, especially for 
amphiphilicity /3 = 21°. In particular the model predicts 
a much smaller, though still finite, critical dipole-field 
strength, agreeing better with our numerical simulations. 
We conclude that interface deformations strongly affect 
the orientation of a tilted Janus particle adsorbed at a 
fluid-fluid interface. 

Finally, Fig. shows that the Janus particle deforms 
the interface in a dipolar fashion: a symmetric depression 
and rise on opposite sides of the particle, as also shown 
in Fig. |2b[ Since the strength of the capillary interac¬ 
tions between two particles interacting as polar capillary 
dipoles depends on the maximal interface deformation 
height ^ 1 , we can tune the strength of capillary interac¬ 
tions by varying the dipole-field strength. These unique 
dipolar capillary interactions may be used to assemble 
particles into novel materials at a fluid-fluid interface in 
a tunable way [TS] . 

V. CONCLUSION 

We studied the behaviour of a magnetic spherical am¬ 
phiphilic Janus particle adsorbed at a fluid-fluid interface 
under the influence of an external magnetic field directed 
parallel to the interface. 

The interaction of the particle with the external field 
results in a torque that tilts the particle, introducing in¬ 
terface deformation. We derived analytical models as¬ 
suming no interface deformations and small interface de¬ 
formations that enable the calculation of the free energy, 
AE, and hence the equilibrium orientation tp, of the par¬ 
ticle in terms of the particle size R, particle amphiphilic¬ 
ity /3, fluid-fluid interface tension 712 , maximal interface 
height Cl, and magnetic dipole-field strength B. We used 
lattice Boltzmann simulations that incorporate interface 
deformations fully to test the results of our analytical 
models. 

In the absence of a magnetic held, our simulations 
showed that the maximal interface deformation height Ci 
increases linearly Ci ~ R'4’ for small tilt angles and then 
reaches a plateau for large tilt angles Ci ~ const. They 
also showed that the free energy of the particle increases 
quadratically for small tilt angles and linearly for large 
tilt angles. We explain this behaviour in terms of our 
theoretical model assuming small interface deformations 
that correctly predicts this behaviour, and hnd that the 


free energy varies as AE = for small tilt angles 

and as AE = 4 i?^ 7 i 2 sin I3ip + C for large tilt angles. 

With the magnetic held switched on, we calculated 
the dependence of the tilt angle ip, on the dipole-held 
strength B, by minimising the calculated free energies 
with respect to the tilt angle tp, for a given dipole-held 
strength B. 

Our undeformed interface model based on free energy 
differences agrees qualitatively with our simulation re¬ 
sults only for large particle tilt angles, but deviates sig- 
nihcantly quantitatively. In particular, this model pre¬ 
dicts a critical dipole-held strength Be = 4 i?^ 7 i 2 sin pi in 
order to rotate the particle, which we do not observe in 
our simulations. Our results contrast with other free en¬ 
ergy difference based models that assume a hat interface, 
which were able to predict the equilibrium orientation of 
nanoparticles reasonably accurately nnun]. 

Our analytical model that includes small interface de¬ 
formations captures the qualitative behaviour of the par¬ 
ticle well for all tilt angles, and also performs well quan¬ 
titatively, in particular for small particle amphiphilicities 
pi. Previous studies by Davies et al. [TS] have shown that 
interface deformations can signihcantly alter the quan¬ 
titative behaviour of ellipsoidal particles at interfaces, 
however, our results demonstrate that interface deforma¬ 
tions can both quantitatively and qualitatively affect the 
behaviour of Janus particles at huid-huid interfaces and 
therefore cannot be neglected. 

Finally, we show that the interface deformations 
around a spherical magnetic titled Janus particle inhu- 
enced by an external held directed parallel to the inter¬ 
face are dipolar in nature. Further, we demonstrate that 
the magnitude of these deformations can be dynamically 
tuned, and therefore the capillary interactions between 
monolayers of such particles are tunable, and we suggest 
that this may enable the production of novel, reconhg- 
urable materials. 


ACKNOWLEDGMENTS 

Financial support is acknowledged from NWO/STW 
(Vidi grant 10787 of J. Harting and STW project 13291). 
We thank the Julich Supercomputing Centre and the 
High Performance Computing Center Stuttgart for the 
technical support and the CPU time which was allo¬ 
cated within grants of the Gauss Center for Supercom¬ 
puting and the “Partnership for Advanced Computing” 
(PRACE). GBD acknowledges Fujitsu and EPSRC for 
funding. 


[1] P. G. de Gennes. Rev. Mod. Phys., 64:645, 1992. 

[2] A. Walther and A. H. E. Muller. Chem. Rev., 113:5194, 
2013. 

[3] J. Yan, M. Bloom, S. G. Bae, E. Luijten, and S. Granick. 


Nature, 491:578, 2012. 

[4] S. K. Smoukov, S. Gangwal, M. Marquez, and O. D. 
Velev. Soft Matter, 5:1285, 2009. 

[5] B. Ren, A. Ruditskiy, J. H. Song, and I. Kretzschmar. 



9 


Langmuir, 28:1149, 2012. 

[6] A. Ruditskiy, B. Ren, and I. Kretzschmar. Soft Matter, 
9:9174, 2013. 

[7] B. P. Binks and P. D. I. Fletcher. Langmuir, 17:4708, 
2001 . 

[8] T. Ondarguhu, P. Fabre, E. Raphael, and M. Veyssie. 
Journal de Physique, 51:1527, 1990. 

[9] D. J. Adams, S. Adams, J. Melrose, and A. C. Weaver. 
Colloid Surface A, 317:360, 2008. 

[10] B. J. Park and D. Lee. ACS Nano, 6:782, 2012. 

[11] B. J. Park and D. Lee. Soft Matter, 8:7690, 2012. 

[12] H. Rezvantalab and S. Shojaei-Zadeh. Soft Matter, 
9:3640, 2013. 

[13] D. Stamou, C. Duschl, and D. Johannsmann. Phys. Rev. 

E, 62:5263, 2000. 

[14] Dominic Vella and L. Mahadevan. Am. J. Phys., 73:817, 
2005. 

[15] J. C. Londet, A. M. Alsayed, J. Zhang, and A. G. Yodh. 
Phys. Rev. Lett., 94:018301, 2005. 

[16] P. J. Yunker, T. Still, M. A. Lohr, and A. G. Yodh. Na¬ 
ture, 476:308, 2011. 

[17] T. Brngarolas, B. J. Park, M. H. Lee, and D. Lee. Adv. 
Funet. Mater., 21:3924, 2011. 

[18] G. B. Davies, T. Kruger, P. V. Coveney, J. Harting, and 

F. Bresme. Soft Matter, 10:6742, 2014. 


[19] G. B. Davies, T. Kruger, P. V. Goveney, J. Harting, and 

F. Bresme. Adv. Mater., 26:6715, 2014. 

[20] S. Sncci. The Lattice Boltzmann Equation for Fluid Dy- 
namies and Beyond. Oxford University Press, 2001. 

[21] X. Shan and H. Chen. Phys. Rev. E, 47:1815, 1993. 

[22] S. Cappelli, Q. Xie, J. Harting, A. M. Jong, and M. W. J. 
Prins. New Biotechnology, in press, 2015. 

[23] A. J. C. Ladd and R. Verberg. J. Stat. Phys., 104:1191, 
2001 . 

[24] F. Jansen and J. Harting. Phys. Rev. E, 83:046707, 2011. 

[25] F. Gunther, F. Janoschek, S. Frijters, and J. Harting. 
Comput. Fluids, 80:184, 2013. 

[26] Y. H. Qian, D. D’Humieres, and P. Lallemand. Euro¬ 
physics Lett., 17:479, 1992. 

[27] S. Frijters, F. Gunther, and J. Harting. Soft Matter, 
8:6542, 2012. 

[28] F. Gunther, S. Frijters, and J. Harting. Soft Matter, 
10:4977, 2014. 

[29] S. Frijters, F. Gunther, and J. Harting. Phys. Rev. E, 
90:042307, 2014. 

[30] C. K. Aidun, Y. Lu, and E.-J. Ding. J. Fluid Mech., 
373:287, 1998. 

[31] S. Dasgupta, M. Katava, M. Faraj, T. Auth, and 

G. Gompper. Langmuir, 30:11873, 2014. 



